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We study the decays B° — > J/i/)7r"'"7r~ and B'^ — > J/ip-K'^n°, including intermediate resonances, 
using a sample of 382 million BB pairs recorded by the BABAR detector at the PEP-II e'^e~ B 
factory. We measure the branching fractions ->■ J/ip p") = (2.7 ± 0.3 ± 0.17) x 10"^ and 

B{B+ J/^ p+) = (5.0 ± 0.7 ± 0.31) X 10"^ 
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We also set the following upper limits at the 90% confidence level: B{B^ — > J/ipir'^'K non- 
resonant) < 1.2 X 10"^, B{B° J/V'/2) < 4.6 X 10"^, and B{B+ J/^tv+tv° non-resonant) < 
4.4 X 10~^. We measure the charge asymmetry in charged B decays to J/ip p to be —0.11 ±0.12 ±0.08. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



The decay B° Q can in principle be used 

to measure the CP violation parameter sin2/3. However, 
the measurement is not as straightforward as for J/i/jK^ 
[3i 9 J because it involves the decay of a pseudoscalar me- 
son to two vector mesons, resulting in both CP-odd and 
CP-even final states. Furthermore, the decay can pro- 
ceed through either a color-suppressed tree diagram, or 
a penguin diagram, both shown in Fig. [l] and interfer- 
ence between them could result in direct CP violation 
Direct CP violation may also occur in 5+ J/4'P~^ de- 
cays, where it would manifest itself as a non-zero charge 
asymmetry: 



A, 



CP 



N{B- J/i^p-) - N{B+ JMjp+) 
N{B- J/tpp-) + N{B+ JMjp+) ' 



(1) 



The large intrinsic width of the p meson necessitates an 
analysis of a significant portion of the invariant mass 
spectrum of the dipion system. 




FIG. 1: Tree and penguin diagrams for the process _B° 



The branching fraction for B''^ —f J/i/'tt+tt^ has previ- 
ously been measured at BaBAR to be (4.6 ± 0.7 ± 0.6) x 
10~^ including a J/ijjp^ component with a branching 
fraction of (1.6 ± 0.6 ± 0.4) x 10^^. This measurement 
used a data sample containing approximately 56 million 
BB pairs, which is a subset of the sample used in this 
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analysis. The charged B decay to J/ipp'^ has not pre- 
viously been observed, the CLEO collaboration set an 
upper limit B{B+ Jj^p-^) < 7.7 x 10~^ at the 90% 
confidence level @. 

The data sample used here contains 382 million BB 
pairs collected with the BaBAR detector at the PEP-II 
asymmetric-energy e~^e~ storage ring, taken at a center- 
of-mass (CM) energy equivalent to the mass of the T{AS) 
resonance. An additional data sample, corresponding to 
an integrated luminosity of 36.8 fb~^, taken at a CM en- 
ergy 40 MeV below the T{4S) resonance, is used to study 
backgrounds from continuum qq production, where q — 
u, d, s, c. 

A detailed description of the BABAR detector can be 
found elsewhere Charged-particle trajectories are 

measured by a five-layer silicon vertex tracker (SVT) 
and a 40-layer drift chamber (DCH) operating in a 
1.5 T solenoidal magnetic field. A detector of internally 
reflected Cherenkov light (DIRC) is used for charged 
hadron identification. Surrounding this is a CsI(Tl) elec- 
tromagnetic calorimeter (EMC), and finally the instru- 
mented flux return (IFR) of the solenoid, which consists 
of layers of iron interspersed with resistive plate cham- 
bers or limited streamer tubes. 

The J/'0 meson is reconstructed in decays to 
where refers to a charged lepton, or /i*. Elec- 
trons are selected on the basis of the ratio of EMC 
shower energy to track momentum, and the energy pro- 
file of the EMC shower. For J/tp — > e~^e~, an attempt 
is made to recover energy losses from bremsstrahlung, 
by looking for showers in the EMC close to those from 
the electron candidates. This procedure increases the 
selection efficiency for J/^p — > e~^e~ candidates by ap- 
proximately 30% ,8] . The muon selection algorithm uses 
a neural network, for which the most important input 
is the number of interaction lengths traversed in the 
IFR. The lepton pairs are fitted to a common vertex 
and the invariant mass of the combination is required to 
be in the range 2.98 (3.06) to 3.14GeV/c2 for the e+e" 
(/i"'"/x~) channels. In order to reduce the background 
from B° — > J/tp K*'^ {K*'^ — > K^Tr~') decays, charged pion 
candidates are required to satisfy stringent particle iden- 
tification criteria, based on combined ionization energy 
loss (dE/dx) in the DCH and SVT with the Cherenkov 
angle measured in the DIRC. 

All tracks are required to originate close to the interac- 
tion point, and to lie in polar angle ranges where particle 
identification efficiency is well measured. The allowed 
ranges correspond to the geometric acceptances of the 
DIRC for pions, the EMC for electrons, and the IFR for 
muons. 

Neutral pion candidates are formed by combining pairs 
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of isolated showers in the EMC. These are required to 
spread over a minimum of three crystals, and to have an 
energy greater than 200 MeV. 

To form a B candidate, the reconstructed J/tp is com- 
bined with cither a pair of oppositely charged pions, or 
a charged pion and a tt^, and a kinematic and geomet- 
ric fit is used to ensure that all final state particles are 
consistent with coming from the same decay point. In 
this fit, we constrain the invariant mass of the and 
the 77 to have the nominal mass of the J/ip and 7r°, 
respectively @. The energy difference, AE, between 
the candidate energy and the single beam energy, , 
(both in the CM frame) is expected to be close to zero 
for signal events, and is therefore required to be in the 
interval -40 to 40 MeV (-60 to 80 MeV) for 5° {B+) 
candidates, corresponding to approximately ±3cr of the 
AE resolution. Note that the range is asymmetric for 5+ 
candidates because the 7r° in the final state gives rise to a 
tail on the low side of the distribution, due to the EMC 
response to photons. For events where more than one 
B candidate passes the selection criteria, the candidate 
with the smallest value of |AE| is chosen. 

The branching fraction for each signal channel is ob- 
tained from: 



B 



X BiJ/ip 1+1-) 



(2) 



where Nsig and Egig are the observed yield and selection 
efficiency, respectively, for a specific signal channel, and 
N is the number of B meson pairs. We assume that 
the T(4S') decays equally often into neutral and charged 
B meson pairs. The J/'0 — > l^l^ branching fraction is 
taken to be (11.87 ±0.12)% 0. 

We extract the signal yields for the J/ipp^ , Jf-iljiT^iT^ 
non-resonant, and Jj"^ ji channels by performing a fit on 
the sample of reconstructed B'^ candidates. We also per- 
form a similar fit to the sample of charged B candidates 
in order to obtain the signal yields for the decay chan- 
nels Jj'^l^P^ and B^ — > J/V'tt+tt^ non-resonant. 
The fits are two-dimensional, extended, unbinned max- 
imum likelihood fits to the distributions of toes and 
TOtttt, Seven event categories are considered: (i) J/V'P sig- 
nal, (ii) J/t/itttt non- resonant signal, (iii) J/V'/2 signal, 
(iv) J/i/' events, (v) background events that do not 
contain a J/V' (non-J/'i/) background), (vi) background 
events containing a J/i/j (inclusive J/i/' background) , and 
(vii) selected background channels that have been stud- 
ied in more detail (exclusive J/'0 backgrounds). In the 
fit to neutral B candidates, the decay channels that com- 
prise category (vii) are J/i^K*^, Jj^K*"-, J/^P Ki{1270), 
J/'4jK+, J/-0P+ [ifll, and J/V'7r+. For the fit to charged 
B candidates, the exclusive J/tp background channels are 
J/^pK*°, J/ipK*+, J/i}Ki{1270), J/i)K+, J/ipK% and 
J/ipK'^. In both cases, these decay channels are not in- 
cluded in category (vi). Of course, categories (iii) and 
(iv) are only present in the fit to neutral B candidates. 

A probability density function (PDF) is constructed 
for each category, and the sum of these PDFs is used to 



fit the data. The likelihood function for the total sam- 
ple is the product of the PDF values for each candidate, 
multiplied by a Poisson factor: 



1 ^ 
£ = ^e-^'(A^')^n^- 



(3) 



where N and N' are the numbers of observed and ex- 
pected events, respectively, and Vi is the value of the 
total PDF for event i. For all event categories except 
for the exclusive J/^' background. Pi is a product of one- 
dimensional PDFs in Toes and TOttu-- 

Fig. [2] shows the toes and to^^t distributions for the 
data, and the projections of the PDFs for each cate- 
gory. The functional forms of these PDFs are as follows. 
For the J/ipp^, J/-07r+7r", J/'!/'/2, and J/tp K'^ compo- 
nents, the Toes distributions are parametrized by Gaus- 
sian functions, all with the same values for the mean and 
width, which are allowed to float in the fit. In the fit 
to charged B candidates, a Crystal Ball function [ll[ is 
used instead for the TOes distributions of the J/tp p+ and 
J/ipn~^TT° signal components, as the presence of a 7r° in 
the final state gives rise to a tail on the low mass side of 
the peak. 

The J/tp p signal component is modeled by a relativis- 
tic P-wave Breit-Wigner function [l2| in TOtttt: 



f + 1 



{{ml 



to2^)2 +TO2r(TO^^)2)' 



(4) 



The pa- 



where r(TO^^) = ) \T+W^ ^ 

rameter q{mTrTr) is the pion momentum in the dipion rest 
frame, with qo = q{mp); P is the J/ip momentum in the 
B rest frame; Le/ / is the orbital angular momentum be- 
tween the J/ip and the p which can be 0, 1 or 2; i? is 
the radius of the Blatt- Weisskopf barrier factor [ij, [l3| , 
which is taken to be (0.5 ±0.5) fm, and nip is the p meson 
mass. 

The TOtttt distribution for the J/^pmr non-resonant sig- 
nal is Ftttt = q{m^Tr)P^ , the product of a three-body 
phase space factor q{m.^Tr)P and a factor motivated 
by angular momentum conservation. 

For the J/V'/2 component, the TOtt^t distribution is de- 
scribed by a relativistic D-wave Breit-Wigner, similar to 
Eq. m but with an extra factor {q/qof' in the expression 
for r(TO^^). 

The decays to J/ip are not considered signal for this 
analysis. Most of them are removed by the requirement 
that all tracks are consistent with coming from the same 
vertex. The to^^t distribution of the remaining J/jpK^ 
events are modeled by a narrow Gaussian function. 

Non-J/ip back grou nd events are modeled by an 
ARGUS function [1^ in toes. The to^^ PDF is the sum 
of two Weibull functions [la] , and a Breit-Wigner to de- 
scribe the p component of the continuum background. 
The parameters of this PDF are fixed to values obtained 
from fits to the J/4> mass sidebands of the data. 



6 



The mES distribution of the inclusive J/ip background 
is an ARGUS function plus a Gaussian at the B mass. 
The width of this Gaussian is somewhat wider than that 
used for signal components as it represents B candidates 
that are not correctly reconstructed. The m^^ PDF is a 
4th-order polynomial. The PDF parameters for this com- 
ponent are fixed to values obtained by fits to a large sam- 
ple oiB ^ J/^{-^ l+l-)X Monte Carlo (MC) simulated 
events, with signal events and exclusive J/tp background 
channels removed. 

Each of the exclusive J/ip background channels is mod- 
eled by a two-dimensional PDF derived from the distri- 
bution of MC events for that decay channel. The nor- 
malizations of these PDFs are determined by taking into 
account the selection efficiency on MC simulation, and 
the world average branching fractions 9] . 

For the branching fraction fit to neutral B candidates 
there are twelve free parameters: the yields of the J/i]^ 
J/ijjTnr, J/ipf2, J/tpK'^, and inclusive J/tp background 
components, the mean and width of the Gaussian used 
for the signal distribution in toes, the parameters nip, 
Fq, and Leff in the p lineshape, and the mean and width 
of the TO^TT distribution for the J/ijjK^ component. All 
other parameters are fixed, including those describing 
lineshape of the /2(1270), the normalizations and shapes 
of the exclusive J/tp background PDFs, and the shapes 
of the inclusive J/tp and non- J/ip background PDFs. We 
also fix the ratio of non- J/ip to J/ip (inclusive plus exclu- 
sive) background yields to a value obtained from fitting 
to data in the region toes < 5.26 (i.e. lower in mass than 
the signal region), and extrapolated to the fit region us- 
ing distributions from MC simulation. 

The configuration for the charged B branching fraction 
fit is very similar. Here, there are eight free parameters, 
since there are no J/%pK^ or J/tp f2 components. 

We find from MC simulation studies that correlations 
between toes a-nd TOtttt give rise to small biases in the 
numbers of Jf-p p^ and J/'p tt+tt'^ non-resonant signal can- 
didates found in the charged B fit. The sizes of these bi- 
ases are evaluated by examining the distribution of resid- 
uals {Nobs — Ninput) for a large number of MC exper- 
iments, and are listed in Table H] The yields obtained 
from the branching fraction fit are therefore corrected to 
take account of this by subtracting these quantities from 
the fitted yields. 

The signal yields and statistical errors obtained from 
the branching fraction fits are fisted in Table HI We also 
lis t the statistical sig nificances of the observed signals, 
a/- 2 Iii{CnuII I Cm ax ) , where CMax is the likefihood from 
the fit, and CnuU is the value of the likelihood func- 
tion when the fit is performed with the signal yield con- 
strained to zero events. 

We obtain signal efficiencies using samples of MC sig- 
nal events, produced in monthly blocks so as to match 
variations in detector and background conditions. Par- 
ticle identification efficiency is corrected using data con- 
trol samples of electrons, muons, and pions. The sizes 
of these corrections vary with momentum and polar an- 
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FIG. 2: Distributions of (a) mES and (b) m,rir for 
J/ip-K^-K~ candidates. The solid line represents the total PDF, 
while the other lines represent (cumulatively, from the bot- 
tom of the plot) non-J/tp background, inclusive J/ip back- 
ground, exclusive J/tp background, J/i/itt^tt" non-resonant 
signal, and J/V'/2 signal. The points with error bars repre- 
sent the data and statistical errors. Plots (c) and (d) show 
the same distributions for B'^ —+ J/tpTv'^-K^ candidates. The 
sharp spike in (b) corresponds to J/ip events, while the 
broader peak is due to J/ip events. 



gle, and average corrections are about 1.5% for electrons, 
5.9% for muons, and 1.8% for pions. With these correc- 
tions applied, about 85% (50%) of electron (muon) pairs, 
and about 85% of pions, satisfy their respective particle 
identification requirements. A small, energy-dependent 
correction (typically about —2% relative) is also applied 
to decay modes containing a tt^ to account for known 
differences in photon detection efficiency between data 
and MC simulation. The corrected signal efficiencies are 
listed in Table HI 

Systematic errors on the branching fraction measure- 
ments arise from uncertainties on the signal efficiency, on 
the fitted yield, on the number of B^B^ or B^ B~ events 
in the sample, and on the J/ip l^l^ branching fraction. 
The number of BB pairs is known to 1.1% accuracy, and 
an additional 1.6% uncertainty is assigned correspond- 
ing to the assumption that the T{AS) decays 50% of the 
time into B^B^ and 50% of the time into B+B" [§1 The 
fractional uncertainty on B{J/ip is 1.0% Q. 

The systematic uncertainties on the efficiency are 
largely due to imperfect simulation of the detector per- 
formance. These effects are studied using various data 
control samples. The largest sources of uncertainty are 
pion identification efficiency, a 2.0% (3.4%) relative error 
for charged (neutral) B decay channels, and tt" efficiency 
(3% for charged B decays). Tracking efficiency (1.5%) 
and lepton identification efficiency (1.0%) also contribute 
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TABLE I: Signal yields, detection efBciencies, and branching fractions for the signal decay channels. The fit bias, product of 
secondary branching fractions {B{J/ip I'^l') and S(7r° —> 77)), and significances of the signals (using statistical uncertainties 
only) are also listed. The corrected yields are obtained by subtracting the fit bias from the fitted yields. For the yields, 
efBciencies, and branching fractions, the first errors are statistical and the second are systematic. For decay channels where no 
significant signal is observed, we quote an upper limit at the 90% confidence level. 



Mode 


Fit bias (events) Corrected yield (events) 


6(%) 




Signif. 


(a) S(xfO"") 







25f.l ± 27.5 ± If. 2 


20.6 ±0.f ±0.8 


ff.87±0.f2 


f3.0 


2.7 ±0.3 ±0.2 


J/lp TT^TT 





64.5 ±35.5 ±7.7 


20.3 ± O.f ± 0.8 


ff.87±0.f2 


2.0 


< f .2 (90% C.L.) 







24.4 ± f3.8 ± f.8 


20.3 ± O.f ± 0.8 


ff.87±0.f2 


2.0 


< 0.46 (90% C.L.) 




-6.8± f.f 


2f8.5± 28.8 ±9.5 


9.7 ± O.f ±0.4 


ff.73±0.f2 


ff.6 


5.0 ±0.7 ±0.3 




+8.2 ±0.9 


-f2.7±27.f ±4.7 


ff.9±0.f ±0.5 


ff.73±0.f2 





< 0.44 (90% C.L.) 



to the uncertainty on the efficiency. The polarization of 
the p in i? — > J/ip p decays is unknown. We use an MC 
sample in which the p mesons are unpolarized to obtain 
the central value of the signal efficiency. We also evaluate 
the efficiency using MC data samples with different p po- 
larizations, and observe a relative variation of 2%, which 
is assigned as a systematic uncertainty on the branching 
fraction measurement. 

We evaluate the impact of the fit procedure by observ- 
ing the changes in the yields when varying the PDF pa- 
rameters that were fixed in the fit within their uncertain- 
ties. The resulting differences are added quadratically for 
sets of parameters that are relatively uncorrelated, and 
added linearly for highly correlated sets of parameters. 
We also repeat the fit using alternative functional forms 
for some PDFs, namely the shape of the inclusive J/ip 
background in m^7r, and the p lineshape, and include the 
resulting differences in the yield in the systematic uncer- 
tainty. In addition, for the J/V'p^ and J/ipn'^Tr'^ channels, 
systematic uncertainties equal to half of the bias correc- 
tions listed in Table U are assigned. The total systematic 
uncertainties on the yield vary from 1.8 events for the 
J/i^ f2 channel, to 11.2 events for the J/ipp^ channel. 

In order to assess the charge asymmetry A'' , we per- 
form a second fit to the charged B candidate sample. In 
this fit, all the shape parameters for the signal and back- 
ground components are fixed to values obtained from the 
branching fraction fit. This reduces the number of free 
parameters and improves the reliability of the fit. We in- 
clude terms for the asymmetries in signal and background 
components as follows: 



1 



^ NP X -(1 
2^ 

+ iV^«xi(l-Q,^^«)7'f« 



bkg 



(5) 



where iV, N'^^, and Nj'^^ are the yields for the J/t/j p^ 
signal, the J/i/jTr+Tr^ non-resonant signal, and the dif- 
ferent background components j, respectively, Qi is the 
charge of the B candidate in event i, and A^, A^^, and 
A'^'^^ are the corresponding charge asymmetries. The 
asymmetry parameters for the exclusive J/^p background 



channels are fixed to world average values P] . The asym- 
metries for the non-J/'ip background and inclusive J/ip 
background components are assumed to be the same 
{A'^'^9 ^ ^bkg^ ^ ^bkgy rpj^.g therefore has six free 

parameters: the yields of the J/ip p'^ signal, J/ipn'^n'^ 
non-resonant signal, and inclusive J/il) background com- 
ponents, and the asymmetries A^ ^ A^^^ and A!^''^ . 

From the charge asymmetry fit, we obtain = 
—0.11 ± 0.12(stat.). The signal and background yields 
obtained from this fit are entirely consistent with those 
from the branching fraction fit. 

A potential contribution to the systematic uncertainty 
on the charge asymmetry A'' could come from differ- 
ent pion identification efhciencies for 7r+ and 7r~, lead- 
ing to different signal selection efhciencies for positively 
and negatively charged B candidates. Using data control 
samples, this effect is found to be negligible. 

The other sources of systematic error on the asymme- 
try are potential differences in the backgrounds for posi- 
tive and negative B candidates. The parameters describ- 
ing the charge asymmetries of the exclusive J/i/' back- 
ground channels are varied within their uncertainties [9|] , 
assuming a 10% uncertainty for the J/il) Ki{127Q) channel 
for which no measurement is available. The normaliza- 
tions of the exclusive background channels, and the shape 
parameters of the inclusive J/i}) background and non-J/^p 
background components are varied in turn, and the fit is 
repeated. The resulting changes to the fitted value of 
A'' are added in quadrature, and the total systematic 
uncertainty is found to be ±0.08. 

In summary, we measure the following branching frac- 
tions, where the first error in each case is statistical and 
the second is systematic: B{B° J/ipp°) = (2.7 ± 0.3 ± 
0.2) X 10-^ and B{B+ J/iP p+) = (5.0 ± 0.7 ± 0.3) x 
10^^ The signals for S" ^ J/V'/2, J/^jh+tt- 
non-resonant, and _B+ J/ipn^Tr'^ non-resonant are not 
statistically significant, thus we set the following upper 
limits at the 90% confidence level: B{B° J/^P 12) < 
4.6 X 10-^ B{B° J/ip-K+TT-) < 1.2 x 10-5, and 
B{B+ J/V'TT+Tr") < 4.4 X lO"**. These values are cal- 
culated by summing the statistical and systematic un- 
certainties in quadrature, multiplying the result by 1.28, 
and adding it to the central value of the branching frac- 
tion. We measure the charge asymmetry defined in Eq. [1] 
for the decays B± -> J/^p^, Ap = -0.11 ± 0.12 ± 0.08. 
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The ARGUS function can be written as 



A{m) = yy/T^eMa^ - V^)) 
A{m) = 



y < 1 
y > 1 



where ^ is a shape parameter, y = m/mmax and nimax 
is a kinematic limit, equal in this case to half the total 
CM energy. 
[16] The WeibuU function can be written as 

W{m) = CV(m- Mon)^''"'^ exp[-T/(m - Mn.ax)'^], 
where V = {C - l)/(C(Afn,ax - M^nf). 
Mmax is the position of the function maximum, Mon is 
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